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Abstract. In 13 14], we proposed a method to characterize jointly self-similarity 
and anisotropy properties of a large class of self-similar Gaussian random fields. We 
provide here a mathematical analysis of our approach, proving that the sharpest way 
of measuring smoothness is related to these anisotropies and thus to the geometry of 
these fields. 
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1 Introduction and motivations 

In numerous modern applications (geography [10j . biomedical imagery (see [3] 
for example) , geophysics [T3] , art investigation Q] , . . . ) , the data available 
for analysis consist of images of homogeneous textures, that need to be char- 
acterized. For such images, a key issue consists first in describing, within a 
suitable framework, the anisotropy of the texture, and then in defining regu- 
larity anisotropy parameters that can actually and efficiently be measured via 
numerical procedures and further involved into e.g., classification schemes. 

Furthermore, in many cases the analyzed textures display at the same time 
self-similarity and anisotropy properties. This is for example such the case 
in medical imaging (osteoporosis, muscular tissues, mammographies,...), cf. 
e.g. |4l5j . hydrology [12], fracture surfaces analysis [5], 

In [13114) . we proposed a method for studying jointly selfsimilarity and 
anisotropy in images by focusing on a specific classic class of Gaussian anisotropic 
selfsimilar processes. We consider £ p norms of hyperbolic wavelets coefficients 
which permits the use of different dilation factors along the horizontal and ver- 
tical axis. We prove empirically that these £ p norms are maximal for a specific 
ratio of the horizontal and vertical axis, directly related to the anisotropy of 
the model. 

In [2J, we proved that there is a close relationship between £ p norms of the 
hyperbolic wavelet coefficients of a function and its norm in a convenient class of 
anisotropic functional spaces : the anisotropic Besov spaces. The consistence 
properties of the estimators introduced in [13114) can then be reformulated 
using sample paths properties of the model in anisotropic Besov spaces. The 
aim of this paper is to provide mathematical foundations of the empirical results 



of [13114] . Here, we prove what we call a strong optimality result, namely that 
the critical exponent of the model in anisotropic Besov spaces is maximal when 
the parameters of the analyzing space fit these of the analyzed textures, which 
is the exact transcription into a mathematical way of the empirical results 
of plM) . 

The paper is organized as follows. In Section [21 we first present the studied 
self-similar anisotropic model. Thereafter in Section [31 we recall some basic 
facts about anisotropic Besov spaces. In Section 0] we then state our main 
result. The proofs are postponed in Section O 

In what follows, we denote £ + the collection of 2 x 2 matrices whose eigen- 
values have positive real parts. 



2 Presentation of the studied model 

In [13114] , we choosed to investigate the properties of a large class of anisotropic 
Gaussian self-similar fields, introduced in [415] . referred as Operator Scaling 
Gaussian Random Field, in short OSGRF. 

For any matrix Eq belonging in £ + such that Tr(Eo) = 2 and any Hq g 
(0,min\ e sp(E o ) R- e (^))? this class can be defined using the following harmoniz- 
able representation: 

X PtEo ,H (x)= [ (e i ^i>~l)p(C)- (Ho+1) ^(i), (1) 



where x — (xi,x 2 ), £ = (£1,^2) and dW(£) stands for a 2D Wiener measure. 
Here p is an Eq -homogeneous continuous positive function, that is satisfies the 
following homogeneity relationship /o(a B °£) = a/(£) on E 2 . We also assume 
that J(l A |C| 2 )p(C)~ 2 ' ffn+1 '^C < +°° which ensures the existence of the Gaus- 
sian field X (see [5]). 

With this construction, the OSGRF X P: e ,h has stationary increments. 
Further it satisfies an anisotropic version of the scaling property with parameter 

Hq (where = denotes equality for all finite dimensional distributions): 

Va > 0, {X p (a E °x)} = {a H °X p (x)} . (2) 

(as usual a E ° = exp(£ log(a)) = £fc(loga) fc ££ /kl). 

Constructions of i? -homogeneous continuous positive function have been 
proposed in [5] via an integral formula (Theorem 2.11). An alternative con- 
struction, more fitted for numerical simulations, can be found in [5J. 

Let us now give a concrete example of OSRGF. Consider the case where 

£0=1?, " ) with a £ (0, 2) and set 



2-a 



P^i,^) = \^\ L/ao + \^ 



V«o 1 It |l/(2-a ) 



The function p is obviously an E homogeneous positive continuous fonction. 
For any H £ min(a ; 2 — ao), the associated Gaussian field will be 

r e ifei> _ 1 , 

X P ,E ,H (XX,X 2 ) = (|a|1/Qo + | 6 | 1/(2 _ ao)) H + l d ^(0 > 



The scaling property satisfied by this field is then 

Va > 0, {X PtEotHo (a ao x 1 ,a 2 - ao x 2 )} = {a H °X p ^ Ho {x u x 2 )} . 

3 Anisotropic concepts of smoothness 

Our main goal here is to study the sample paths properties of this class of 
Gaussian fields in suitable anisotropic functional spaces. This approach is 
quite natural (see since the studied model is anisotropic. To this end, 
suitable concepts of anisotropic smoothness are needed. The aim of this sec- 
tion is to give some background about the appropriate anisotropic functional 
spaces : Anisotropic Besov spaces. These spaces generalize classical (isotropic) 
Besov spaces and have been studied in parallel with them. Let D a diagonal- 
izable matrix of £ + with eigenvalues Ai, A2 and associated eigenvalues (ei, e 2 )- 
The anisotropic Besov spaces with anisotropy D can be defined as follows (see 
Theorem 5.8 of QI]) : 

Definition 1. Let (p, q) G [1, +00] 2 , s > and /3 6 E. For any / e L P (R 2 ) set 

n/b« „ (R 2 ,-D) = E ( f \\{Ktfm\\i*t- sq/at - 1 \ logwr^dt 

By definition 

^. g .| log |,(K 2 ,£) = {/ G h p (M. 2 ), \\f\\ 6 . ^ >E) < +^} . 

p,g,\ log l p 

The matrix D is called the anisotropy of the Besov space q , log , p (R 2 , D). 

Remark 1. It is well-known that there exists a strong relationship between 
£' p norms and classical wavelet coefficients (see |16|17j ). In |13|14) . our esti- 
mators are based on hyperbolic wavelet analysis. In [2], we proved that i v 
norms of hyperbolic wavelet coefficients of a functions are related to its norms 
in anisotropic Besov spaces. It is the reason why the natural mathematical 
framework to relate anisotropy and self-similar properties of the model to its 
sample paths properties is this of anisotropic Besov spaces. 

Remark 2. Let D a diagonalizable matrix of £ + . For any a > 0, XD is also a 
diagonalizable matrix of £ + with eigenvalues aAi,aA2 and same eigenvectors 
as D. Hence for any s > 0, B™ q (R 2 ,\D) = B s p q (R 2 ,D). 

Using Remark [2j we deduce that without loss of generality, we can assume in 
the sequel that Tr(D) = 2. We then define 

£t = {Ee£ + , Tr(E) = 2} . 

As it is the case for isotropic spaces, anisotropic Holder spaces C S (M. 2 ,E) can 
be defined as particular anisotropic Besov spaces. 




Definition 2. Let s > 0, f3 G K. The anisotropic Holder spaces C| S log , p (R 2 , D) 
are denned by 

C flog|f>( K2 ,£) = B lo,oo,\ log )^) ■ 

Hence, a bounded function / belongs to C* log ^(R 2 , if and only if for any 
r G (0, 1), G S E (\ ■ \ E ) and x G R d , one has 

\f(x + r E 0)-f(x)\<C o r s \\og(r)f , 

for some Co > 0, that is if and only if its restriction f& along any parametric 
curve of the form 

r>04 r E , 

with G Sjfd • |_e) is in the usual Holder space C, s , ^(R) and ||/e||c s (r) 

I § I I log 1 

does not depend on 0. Roughly speaking, the anisotropic "directional" regu- 
larity in any anisotropic "direction" has to be larger than s. In other words, we 
replace straight lines of isotropic setting by curves with parametric equation 
r>0n r E adapted to anisotropic setting. 

To state our optimality results we need a local version of anisotropic Besov 
spaces : 

Definition 3. Let D G £ + be a fixed diagonalizable anisotropy, 1 < p,q < oo, 
(3 e E, < s < oo and / G Lf oc (R 2 ). 

The function / belongs to q |log|/3 ioc (R 2 ,^) if for any <p G V(R d ), the 
function ipf belongs to B£ q , log , p (R 2 ,D). 

The anisotropic local critical exponent in anisotropic Besov spaces „ (R 2 , E) 
of / G if oc (R 2 ) is then defined by 

a f ,i oc (E,p,q) = sup{s, / G B° q loc (R 2 , E)} . 
4 Statement of our main result 

In what follows, we are given Eo G an< ^ Pe an*E -homogeneous continuous 
positive function, H G (0, min Ag 5 p ( f ; )(Re(A)). 

Our results will be based on a comparison between the topology related 
to pe„ involved in the construction of the Gaussian field {X PEgt E ,Ho( x )}xeR d 
defined by equation (fTJ) and this of the analyzing spaces Bp q (R 2 ,D). To be 
able to compare these two topologies, we also assume that D G ■ 

We characterize in some sense an anisotropy Eq and an Hurst index of the 
field {X PEo ^H (x)}x£m. d -i which is interesting when analyzing anisotropic self- 
similar textures. Combining the results of this paper and these of [2], our 
approach can thus be turned in an effective algorithm for the estimation of 
the anisotropy of self similar textures (see [13114) ). In the case where q = 
oo, Theorem 1 is the exact mathematical reformulation of equation (9) stated 
in [M], on which are based the definition of the estimators of the anisotropy 
and smoothness of the model. Our main result in then both a mathematical 
justification and an extension of the results empirically proved in [13114] : 



Theorem 1. Let (p,q) G [1, +00] . Then almost surely 



(Eo,P,q) = s ^p{ a x„ E , E() ,H ,ioc{D,p,q), D diagonalizable and D € 
In fact, Theorem [T] contains two main results : 

• The critical exponent of the field {X PEgt H ( x )}xGR d m anisotropic Besov 
space -Bp 9 (R d , -D ) where D is the real diagonalizable part of D (see 
Proposition 4.1 of [8] for a definition) equals the associated Hurst index 
Hq. It has already been proved in [5] in a general d-dimensional context. 

• The diagonalizable real part Do of any anisotropy Eq of the field {X PEo _h ( x )}xeR d 
maximizes this critical exponent among all possible analysis matrices. In 
fact, the "best way" of measuring smoothness of the field {X PEo t H ( x )} X £R d 

is to measure smoothness along the "anisotropic directions" r > H> r D °G, 
related to the genuine geometry of the field. 

5 Proof of Theorem CD 

Proposition 5.2 of [8] directly implies that a.s. 

ax PE(rHo ,ioc(E ,P,q) = H . 
We now prove that if D is a diagonalizable matrix of then a.s. 

a XpBo<Ho ,ioc(D,p,q) < H a . 

Denote Do the real diagonalizable part of Eq (see Proposition 4.1 of [18 ). Ob- 
serve now that, by Theorem 1.1 of [7] and Lemma 5.1 of [8], one has a.s. 

a X PEQ ,h (A P, q) = ax PDo , Ho (D, p, q) , 

where pe q , Pd q are respectively two Do- and Do-homogeneous continuous pos- 
itive functions and {X PEot H ( x )}xeR, d y {X PDa ,H ( x )}xeR d the associated Gaus- 
sian fields by ([1]). One can then investigate the sample paths properties of 
{X PDo: H ( x )}xm d instead of this of {X pEoiHo (x)} xeM d. 

Hence from now, we then assume that E equals its diagonalizable real 
part, namely that Eq = Do. Even if D and Do are both diagonalizable since 
these matrices are not commuting we cannot assume that Do and D are both 
diagonal. 

Let ip G T>(M. d ). As above, one may assume supp(<y3) C K = Dd o (0, 1). 
Denote A? < A!] (resp Ai < A2) the two eigenvalues of Do (resp. D) and e?, e§ 
(resp ei,e2) some associated eigenvectors. We exclude the two cases Aj = A°, 
Ai = A2 corresponding to the cases D = Id, D = Id which can be deduced 
from Theorem 4.1 of [8]. Hence the inequalities between the eigenvalues are 
strict. 

Let us fix p G [1,+co]. By assumption [Ho/Xl] = [-ffo/A"] = 1- The proof 
of Proposition 5.5 of [8] provides us with the following equalities 

H = suptf,/ 1 \\{yX){x + tel) - ^X)(x)\\ p LP t^/^i dt < +to} 

= sup^Jo \\(<pX)(x + te%) - (<pX){x)\\ p LP t-W x *- x <lt < +00} . 1 ' 



Define now the two followings indices 



Q'l 



i(p) = sup{ J 8, / UtpXXx + tet) - (<pX)(x)\\ p L ,t-W>*- 1 dt < +w} , 



(p) =su P {/3, / ||(^)(a: + te 2 ) - (^(a)!^*-*^- 1 ^ < +00} . 



"2 



Thereafter set a(p) = mm(ai(p), a 2 (p)). Assume that for some p G [1, +00] 
one has a.s. a(p) > Ho, that is a.s (pX G Bp,p (M 2 , 13) with a(p) > i?o- 
We first need the following lemma : 

Lemma 1. Only five cases are possible : 

1. 

H ( H H \ . ( a x {p) a 2 (p) 

AY =max lAf'AfJ- mm l^T'^r 

2. 



Ho (Ho Ho\ . f ai (p) a 2 (p)\ H Q 



-Ho Ho 
A?' A° 



-Ho f Ho H \ H ( Ho H 

-Jo = max [~Jp' -Jo) >Jo = mm (^o> "^0 I n,in 



onip) « 2 (p) 
Ai A2 



Ho (Ho H \ (a x (p) aa (p) . 



5. 



«i(p) Q2(p) 
Ai A2 



> 



#0 
A? 



Ho Ho 
A?'A° 



> mm 



Ai A 2 
otiip) a 2 (p) 



Ai 



A 2 



#0 



-Ho 
A^ 



#0 -H 
Af' Al" 



Hp -Ho 
A?' A° 



Proof of Lemma Q3 The only point to prove is that the case 



,ai(p) a 2 (p) H ,H Q H . ,a x {p) a 2 (p) H . ,H a H 

max (— . — — ) = — = niax (TF' — ) > mm (—J^~' ~T—> = — = mm (TcT> To ' 



Ai 



A 2 



Ai 



Ai A 2 



A 2 



A> 



A" A 2 



is impossible. Suppose that this relationship holds. Then one has 



ai{p) _ -Ho and a 2 (p) _ H \ ^ / a 2 (p) _ Ho ^ a\(p) _ Ho 
Ai A® \ 2 A 2 / \ A 2 A® Ai A2 



\i a 1 A 2 a 2 j \ A 2 

Since min(ai(p),a 2 (jj) > Hq, it implies that 



Ho „ -Ho Ho H \ ( H H Hq Hq 
— < -rp and — < -Q- or _ < _ and — < 

Al Ai A2 Ao / \ A2 Ai Ai A 2 



that is 

(Ai > A? and A 2 > A") or (A 2 > A? and A x > A 2 ) . 

Since A? + A 2 = Ai + A 2 = 2, it yields to a contradiction. 

We now deal successively with the five cases of Lemma[TJ Let us first assume 
that point (f) of Lemma [T] holds. 

Let s G (Hq/\1, min(ai(p)/Ai, a 2 (p)/A 2 , 1)). Then, by definition of s one 

has 

/ \\(<pX)(x + tei) - {vX){x)\\ p LP t- s v-\\t < +00 , 
Jo 

and 

/ \\(<pX)(x + te 2 ) - (<pX)(x)\\ p Lf t--*- 1 dt < +00 . 
Jo 

The, using the finite difference characterization of classical Besov spaces, one 
deduces that ipX e B° p (R 2 ). Since, by Proposition 5.6 of 0, B s pp (R 2 ) c — > 

-Bp^ 1 (M 2 , Z?o) and s > ifo/A? it yields to a contradiction. 

Now suppose that point (2) of Lemma[T]holds. Assume that min(ai(p)/Ai, a 2 (p)/A 2 ) = 
ai(p)/Ai. The other case will be similar. Observe that that if ef = ae\ for 
some agM, one then has 

Ho/A? - sup{ s , / \\(<pX)(x+te 1 )-(i P X){x)\\ p LP t- sp - 1 dt < +^} - a x (p)/Ai , 
Jo 

which is impossible since by assumption Hq/\\ > a\{p)/\i. 
Then the family (ei, e°) is necessarily a basis of R 2 . Hence 

e° = ae x + be , , (4) 

for some (a, 6) el 2 . By for any e > one has a.s. 

/ \\(<pX)(x + te° 2 ) - ^X)(x)f LP t-^P/^+^dt = +00 . (5) 
Jo 

Further, the triangular inequality and equation Q implies that 

fi \\(<pX)(x + te°) - (<fiX)(x)\\%t-( H °^°+*)-idt 
< 2P- 1 Si \\(<pX)(x + t(aei + fee?)) - (<pX)(x + tbe\)\\ p Lv t-^/^ + ^- l dt 
+2P- 1 Si \\(ipX)(x + tbe° 1 ) - (ipX)(x)\\P LP t- ( - H «P/ x "+^- 1 dt . 

(6) 

Set now y = x + ibe\ and remark that 

\\{<pX){x + t(aei + be ,)) - (<pX)(x + tbe\)\\ p LP = \\(<pX)(y + taei) - ( V X){y)\\ p LP 
In addition, H /\°, < «i/Ai < H$/\\. Since 

H a /Xl = sup{s, f \\(<pX)(x + tei) - (<fiX)(x)\\ p LP t- s P- 1 dt < +00} , 
Jo 



and by definition of a.\(p), one deduces that for any e > sufficiently small 

\\{<pX){y+tae x )-{vX){y)\\ p LP t ^ dt < \\(<pX)(y+tae{)-(<pXM\\ 

Jo 

(7) 
and 

/ ||(^0(x+t6e°)-(^0(z)ll£^ ^ + dt< / IK^^+^e?)-^)^)!! 
Jo Jo 

are finite. Combining this information with Equations ©, ©, implies a con- 
tradiction. The proof of Theorem [1] in the three other cases of Lemma [1] is 
exactly similar. 
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